Structure and conductance histogram of atomic-sized Au contacts 
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Many experiments have shown that the conductance histograms of metalhc atomic-sized contacts 
exhibit a peak structure, which is characteristic of the corresponding materiaL The origin of these 
peaks stiU remains as an open problem. In order to shed some light on this issue, we present a 
theoretical analysis of the conductance histograms of Au atomic contacts. We have combined clas- 
sical molecular dynamics simulations of the breaking of nanocontacts with conductance calculations 
based on a tight-binding model. This combination gives us access to crucial information such as 
contact geometries, forces, minimum cross-section, total conductance and transmission coefficients 
of the individual conduction channels. The ensemble of our results suggests that the low temper- 
ature Au conductance histograms are a consequence of a subtle interplay between mechanical and 
electrical properties of these nanocontacts. At variance with other suggestions in the literature, our 
results indicate that the peaks in the Au conductance histograms are not a simple consequence of 
conductance quantization or the existence of exceptionally stable radii. We show that the main 
peak in the histogram close to one quantum of conductance is due to the formation of single-atom 
contacts and chains of gold atoms. Moreover, we present a detailed comparison with experimental 
results on Au atomic contacts where the individual channel transmissions have been determined. 

PACS numbers: 73.63.Rt, 73.40.Jn, 05.10.-a, 02.70.Ns 



I. INTRODUCTION 



The metallic nanowires fabricated by means of the 
scanning-tunneling microscope and break-junction tech- 
niques have turned out to be an unique playground where 
to test concepts on electronic transport at the atomic 
scalei. The great activity in this field during the past 
decade was triggered by the observation of a step-Uke 
evolution of the conductance during the formation of 
the nanowires^^. Immediately, the connection with con- 
ductance quantization (CQ) observed in two-dimensional 
electron-gas devices^ was discussed, where every conduc- 
tion channel has a perfect transparency. In order to in- 
vestigate in detail CQ, different authors introduced con- 
ductance histograms, constructed from a large number 
of individual conductance curve a^i^i'^ , i.e. the evolution 
of the conductance with the elongation. These demon- 
strated that for metals like gold and sodium the con- 
ductance has a certain preference to adopt multiples of 
Go = 2e^//i, the quantum of conductance. However, for 
a large variety of metals, the peaks in the histograms 
do not appear at multiples of Go (for a detailed discus- 
sion of the conductance histograms, see section V.D in 
Ref. [l|). For instance, it has been shown^ that an Al 
contact sustains up to three conduction channels in the 
last conductance plateau, although the typical conduc- 
tance is close to 1 Go. The Al histogram exhibits a series 
of peakg^, which in view of the results of Ref. Q can- 
not be interpreted as CQ. These results showed that one 
cannot simply interpret the peaks close to multiples of 
Go in the histograms as evidence for CQ^. The results 
of Ref. 8] were explained on the basis of a tight-binding 
modelifliii. From this work the picture that emerged 
was that the number of conduction channels is limited 



by the valence orbitals of the corresponding metal, and 
the channels are in general not fully transparent. This 
view was confirmed experimentally by a systematic study 
for various metals^. In spite of the success of this atomic 
picture, the origin of the peaks in the conductance his- 
tograms remains unclear, and the main goal of this work 
is to shed some light on this issue. 

It is obvious that the mechanical properties of the 
nanowires can play an important role in the shape of 
a conductance histogram. This was pointed out in early 
molecular-dynamics simulations^'^'^'*'^^, which suggested 
that the steps in the conductance during the formation of 
the nanowires are due to atomic rearrangements. Clear 
experimental evidence of this fact was reported by Rubio 
et alw^, who combined conductance and force measure- 
ments to show that the jumps in the conductance are 
associated with distinct jumps in the force. A striking 
example of the interplay between mechanical and elec- 
trical properties can be found in the work of Yanson et 
al}'^-^^ . These authors observed a series of peaks in the 
histograms of alkali metals at relative high temperatures, 
which were convincingly interpreted as the existence of 
exceptionally stable contact diameters due to electronic 
and atomic shell effects. This confirmed the idea that the 
quantum modes in nanowires not only determine the con- 
ductance but also play an important role in the cohesive 
energ}ii2iSSiSl. More recently, Hasmy et alm^ have cal- 
culated histograms of the minimum cross-section for Al 
contacts using molecular dynamics. At low temperatures 
they obtained peaks at multiples of the cross-section of 
a single atom, which led them to a new interpretation 
of the conductance histogram peaks based on preferen- 
tial geometrical arrangements of nanocontact necks. It 
is also important to mention that very recently the room 
temperature conductance histograms of Au and Al have 



been interpreted as an evidence of electronic and atomic 
shell effect in these metals^iSi. 

So far, several publications have analysed the evolu- 
tion of metallic nanowires in a stretching even1iSSiS&. The 
theoretical analysis of conductance histograms, however, 
has only been carried out within the framework of free- 
electron models by choosing particular nanowire dynam- 
ical. Some authors have combined classical molecular- 
dynamics simulations with atomistic calculations of the 
conductanceiSiSSiSSiffii^iiSSiSiSi to analyze single realiza- 
tions of the contact. The existence of some exceptionally 
stable structures and the correlation of conductance his- 
tograms and atomic structure were previously published 
for gold using the Wulff-construction method and the 
Hiickel approximation for the conductance^^*^. How- 
ever, there are no full atomistic studies of the conduc- 
tance histograms. 

In this work we report for the first time a full atom- 
istic study of the conductance histograms of gold atomic 
contacts, and we present an extensive comparison with 
experimental results. For this purpose, we have com- 
bined atomistic calculations of the conductance within a 
tight-binding approximation with realistic molecular dy- 
namics simulations of contact geometries. This combi- 
nation allows us to obtain detailed information on the 
mechanical and transport properties such as contact ge- 
ometries, forces, minimum cross-section, total conduc- 
tance, the number and the evolution of individual trans- 
mission coefficients. We have chosen the element gold 
for two reasons. First, gold histograms have been widely 
investigatec&Si2L^2i2242iii. One typically observes a pro- 
nounced peak at 1 Gq, and much smaller peaks close 
to 2 Go and 3 Go, although the histograms slightly de- 
pend on details such as temperature, voltage bias and 
vacuum conditions. Second, the experimental analysis of 
the individual conduction channels using the technique 
introduced in Ref. j3| is availableiS^SiM^i. This anal- 
ysis shows that conduction channels disappear one by 
one as the contact is elongated, and in the last plateau, 
which corresponds to either a single-atom contact or 
a chain^^*^, the conductance is dominated by a single 
channel. 

Our results show that the histograms of the mini- 
mum cross-section at low temperatures (4 K) exhibit 
peaks, suggesting the existence of some exceptionally sta- 
ble structures. From a semi-classical point of view, one 
would then expect these peaks to be reflected as peaks 
in the conductance histograms. However, our results in- 
dicate that this is in general not the case. Indeed, only 
in the last stage of the stretching process we obtain a 
large peak close to 1 Go, which is correlated to the for- 
mation of a single-atom contact or a chain of atoms. Our 
channel analysis shows that the conduction modes disap- 
pear one by one and the last plateau is dominated by a 
single channel, in agreement with the experimental re- 
sults. In short, the ensemble of our results shows that 
the low-temperature conductance histograms of gold are 
a consequence of the interplay between the mechanical 



and the electrical properties of these atomic contacts. 

The rest of the paper is organized as follows. In sec- 
tion II we explain the technical details of our molecular 
dynamics simulations and conductance calculations. Sec- 
tion III is devoted to the analysis of some examples of the 
breaking of the contacts, where in particular we show the 
evolution of the single channels and illustrate the typical 
geometries formed in the last stages: single-atom con- 
tacts and chains of atoms. In section IV we present new 
experimental results of the determination of the channel 
transmissions in Au contacts and we compare these re- 
sults with our theoretical simulations. In section V we 
discuss our results for the histograms of the minimum 
cross-section. The corresponding results for the conduc- 
tance histogram are presented in section VI and their 
correlation with the minimum cross-section histograms 
is discussed in detail. Finally, in section VII we present 
a brief summary of the main results and conclusions. 



II. DESCRIPTION OF THE THEORETICAL 
MODEL 

As mentioned in the introduction, in order to study the 
histograms of gold nanowires we have performed classi- 
cal molecular dynamics (MD) simulations to investigate 
the formation of the contacts, and we have used a tight- 
binding model to calculate the conductance during the 
breaking of the atomic necks. It is worth stressing that 
the resulting numerical effort is considerable . In what 
follows we explain the details of both theoretical tech- 
niques. 



A. Molecular dynamics simulations 

We analyze the breaking of Au nanocontacts by means 
of classical MD simulations. The forces and energies are 
calculated using semi-empirical potentials derived from 
the effective medium theory^*, which successfully de- 
scribe experimental bulk and surface properties'*^ and 
have also been successfully used for simulating nanowires 
by Jacobsen et. al.2ii^2i2ia2&. 

Transmission electron microscopy^ shows that in the 
last stage of the stretching process, nanocontacts of Au 
are crystalline and atom rearrangements take place in 
such a way that one of the crystal directions [100], [110] 
or [111] lies in stretching direction of the wire, indepen- 
dent of the starting crystal orientation. For that reason 
our starting configuration is a perfect fcc-lattice of 112 
atoms in [100] direction (z direction) of length 2.65 nm 
with a cross-section of eight atoms. This central wire 
is attached at both ends to two slabs of 288 atoms each, 
which are kept fixed. We also performed simulations with 
starting configurations of perfect lattices in [110] (central 
wire: 108 atoms; each slab: 196 atoms) and in [111] di- 
rection (central wire: 240 atoms; each slab: 720 atoms). 
The stretching process is simulated by separating both 



slabs symmetrically at fixed distance in every time step. 
The resulting stretching velocity of 2 m/s is much bigger 
than in the experiments, but it is small compared with 
the speed of sound in the material (3000 m/s). Thus 
the wire can re-equilibrate between successive instabil- 
ities, although collective relaxation processes might be 
suppressedi^i^. The Newtonian equations of motion of 
the wire atoms are integrated via the velocity Verlet al- 
gorithm^ with a time step of 1.4 fs. A Nose-Hoover 
thermostatic maintains the average temperature at 4.2 
K. We used periodic boundary conditions for the slabs in 
z direction and minimum image convention for the slabs 
perpendicular to the z direction^. Before the stretching 
process, every atom of the wire gets a randomly chosen 
velocity (that is why every nanocontact evolves a little 
bit different from stretching to stretching) and the wire 
is equilibrated for about 0.7 ns with periodic boundary 
conditions perpendicular to the z direction. 

In order to test whether the conductance changes are 
correlated with atom rearrangements in the nanocontact, 
we calculate the radius of the minimum cross-section 
perpendicular to the stretching direction as defined in 
Ref. |23|. For this purpose, every atom is represented by 
a sphere with the volume of the elementary cell in the 
fcc-lattice. For a given configuration a slice with a width 
of about the interlayer distance^ is taken perpendicu- 
lar to the stretching direction. Out of the volume of the 
atomic spheres overlapping with the slice, the radius of 
a cylinder which would fill the same volume in that slice 
is computed. The procedure is repeated along the whole 
nanocontact and the smallest radius is taken as the radius 
of the minimum cross-section of a given configuration. 

Moreover, during the stretching process, every 1.4 ps 
a configuration is recorded and the strain force of the 
nanocontact is computed following Finbow et al^. Ev- 
ery 5.6 ps the corresponding conductance is calculated 
using the method described in the next subsection. 
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where i,j run over the atomic sites and a,P denote the 
different atomic orbitals, eia are the on-site energies and 
Via.jp are the hopping elements. Additionally, we need 
the overlaps between the different orbitals, Siajp- We 
take all the se p arameters from the bulk parameteriza- 
tion of Ref. [53 , which is known to accurately reproduce 
the band structure of bulk materials. The atomic basis 
is formed by the 9 orbitals: 5d, 6s, 6p, which give rise to 
the main bands around the Au Fermi energy. It is im- 
portant to emphasize that in this parameterization both, 
the hopping elements and the overlaps, are functions of 
the distance between the atoms, which allows us to ap- 
ply it in combination with the MD simulations. These 
functions have a cut-off radius that encloses up to five 
nearest neighbors in a bulk geometry. 

In an atomic contact the local environment in the neck 
region is very different to that of the bulk material. In 
particular, this fact can lead to large deviations from the 
approximate local charge neutrality that typical metallic 
elements exhibit. We correct this effect imposing local 
charge neutrality in all the atoms of the nanowire through 
a self-consistent variation of the diagonal parameters e^Q • 

The quantity that we are interested in is the low- 
temperature linear conductance which can be expressed 
in terms of the Landauer formula 



G = GoTr{i{EF)i\EF)} = GoJ2 ^»' 
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where t is the transmission matrix of the contact, and the 
T„'s are the transmission eigenvalues at the Fermi energy 
Ep . The transmission matrix can be calculated in terms 
of the Green functions of the contact as follows 
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B. Conductance calculations 

The computation of conductance in atomic systems of 
the size that we consider here, is out of the scope of 
present ab-initio methods. For this reason, we chose to 
use a tight-binding model based on the parameterization 
introduced in Ref. [53. As mentioned in the introduc- 
tion, the tight-binding models have been very success- 
ful in the description of the transport of metallic atomic 
contactsifliii. Our approach follows closely the one of 
Refs. jirillj |. with the only difference that in this case we 
use a non-orthogonal basis, which introduces some minor 
changes as explained in Refs. |56i57ij . 

The first step in our approach is the description of 
the electronic structure of the atomic contacts, which is 
done in terms of the following tight-binding Hamiltonian 



t{E)^2r'/\E)G^ciEKii (E) 
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Here, F^^^ are the scattering rate matrices given by 

Fl,/? = Im(SL.fl), where Sl,/j are the self-energies which 
contain the information of the electronic structure of the 
leads and their coupling to the central part of the con- 
tact. In our case the leads start in the slabs with 288 
atoms and the central part, denoted by C, is the central 
wire. The self-energies can be expressed as 
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where v is the hopping matrix which describes the con- 
nection between the central cluster and the leads, and 
S is the corresponding overlap matrix. In this expres- 
sion, gL,R are the Green functions of the uncoupled leads, 




FIG. 1: (Color online). Formation of a dimer configuration 
of Au in the [100] direction at T = 4.2 K. The upper panel 
show the strain forces as a function of the elongation of the 
contact. In the lower panel we show the conductance (solid 
line), the radius of the minimum cross-section (dashed line), 
and the channel transmissions (lines in the bottom of the 
graph). The vertical lines define regions with different number 
of open channels ranging from 8 to 1. Below the graph we 
show snapshots of the breaking of the contact. 



which we assume to be the bulk Green functions of Au. 
Finally, the central cluster Green functions, Gc, are given 
by 



Gc{E)= ESc-Hc-tL{E)-tR{E) 



(5) 



where He and Sc are the Hamiltonian and the overlap 
matrix of the central part, respectively. 



III. EXAMPLES OF THE CONTACT 
FORMATION 

In this section we discuss in detail some of the typical 
features that we observe in single realizations of the Au 
contacts. In particular, we discuss the formation of both 
single-atom contacts and chains. 



A. Single-atom contacts 

In Fig. n we show the formation of what we could call 
a single-atom contact. To give a complete information 



of this process we have depicted the forces, the conduc- 
tance, the radius of the minimum cross-section, and the 
transmissions of the different conduction channels. In 
this example the conductance evolves during the elonga- 
tion process from a value close to 5 Gq for the starting 
point to 1 Go in the last stages before breaking^. The 
decrease of the conductance follows closely the evolution 
of the minimum cross-section. This correlation is par- 
ticularly pronounced in the last stages when the contact 
has one or two atoms in the narrowest part. However, it 
is worth stressing that the conductance does not always 
follow the minimum cross-section, as can be seen in the 
elongation region between 0.75 and 1.0 nm. Sometimes 
there appear jumps in the conductance in regions where 
the minimum cross-section evolves smoothly. We inter- 
pret this fact as rearrangements of atoms away from the 
narrowest part of the wire. This interpretation is con- 
firmed by the analysis of the forces, as we explain in the 
next paragraph. 

As shown in the upper panel of Fig. ^ one can see in 
the forces a series of regions where they increase linearly 
(elastic stages) and regions with abrupt jumps (plastic 
stages). The first ones correspond to situations where 
the structure remains fundamentally unchanged, and the 
latter correspond to the breaking of bonds and the sub- 
sequent sudden atomic rearrangements. As commented 
in the introduction, this evolution of the forces has been 
measured for Au contacts at room temperature with the 
help of an atomic- force microscope^^. The order of mag- 
nitude of the forces in Fig. ^ is in good agreement with 
this experiment. Notice that also in the cases where the 
jumps in the conductance are not correlated to abrupt 
changes in the minimum cross-section, one sees jumps in 
the forces. This indicates that plastic deformations in 
regions away from the narrowest part of the wire have 
also an influence on the conductance. Thus, we see that 
in the determination of the conductance, the minimum 
cross-section is an important ingredient, but it is by no 
means the only one. This fact will be crucial later on for 
the interpretation of the conductance histograms. 

Another interesting feature in Fig.^is the last conduc- 
tance plateau. It corresponds to a contact with the cross- 
section of one atom, but it is formed by a dimer of gold 
atoms. This dimer configuration is the most common 
geometry that we find in the last stages of the contact 
breaking. This type of configuration has been reported 
recently in MD simulations of the breaking of Al wiresS. 
Notice also that in the last plateau, which is marked by 
a linear increase of the forces, the conductance is close to 
1 Go and it is clearly dominated b y a sin gle conduction 
channel. As explained in Refs. |l0lllll2| . this is due to 
the fact that the number of channels is controlled by the 
number of valence orbitals in the narrowest part of the 
contact. In the case of Au the main contribution to the 
density of states at the Fermi energy comes from the 6s 
band, i.e. there is a single valence orbital. This implies 
that in the dimer configuration of Fig. ^ there is only 
one possible path (conduction channel) which proceeds 



through the s orbitals. It is also important to emphasize 
that during the elongation of the contact we observe that 
the channels disappear one by one in a gree ment with the 
experimental results reported in Ref. [l^l- This will be 
discussed in more detail below. Let us point out that 
we consider a channel being closed if its transmission is 
below 0.01 Gq. We attribute the successive closing of the 
channels to the fact that in the stretching process the 
Au atoms leave the narrowest part one by one, and every 
Au atom contributes to the transport with one orbital^, 
which in turn can give rise to one conduction channel. 



B. Atomic chains 

In some occasions we observe that the contact does 
not break as in the configuration of Fig. ^ but con- 
tinues to form chains of several atoms. This is illus- 
trated in Fig. 121 Experimental evidence of the forma- 
tion of gold chains was first reported independently by 
two groups^^*^. The formation of chains had been al- 
ready suggested by different simulations^iiSSi^. In the 
last years there has been an intense activity in this topic. 
Thus for instance, there are now several experiments con- 
firming the existence of gold chainsSLSSiS^; chains have 
been also observed in other materiak^"* , and the forces 
during the formation of a Au chain have been mea- 



suredSi. From the theoretical side, many authors have 
analyzed the formation, stability and conductance of gold 

^]^^j^^25.26.57.66.67.68.69.70.71.72.73.74.75.76 p^^ ^ ^^^^ ^^^_ 

plete list of references, see section XI in Ref. [l[. We want 
to point out that so far there is no theoretical analysis of 
the conductance during the chain formation available. 

In our simulations for the [100] direction we have ob- 
served 11 chains out of 50 stretching processes with differ- 
ent lengths: nine chains ranging from 3 atoms to 6 atoms 
and three chains of 10, 12 and 14 atoms. Experimentally 
chains up to 8 Au atoms have been reported^. The 
mechanism of the chain formation has been explained by 
Bahn and Jacobsen2& in terms of many-body effects in 
metals. The main idea is that in certain metallic systems 
the binding energy per neighboring atom may increase 
as the number of neighbors decreases. Here, we want to 
illustrate how a chain is formed in real time. We show in 
Fig. five snapshots of the dynamics of the nanocontact 
of Fig. 121 At 0.0 ps an one-atom contact is shown with a 
typical dimer structure in the middle. The atom number 
2 has four nearest neighbors and the atom 1 only three. 
Due to the many-body interactions the bond between the 
atoms 2 and 3 breaks. After about 40 ps another bond 
breaks and the atom 2 moves into the chain. The atom 4 
is now in a similar situation as the atom 2 was at the be- 
ginning, and the previous process can be repeated again. 

Turning now to the transport properties of the chains, 
we see in Fig. (21 that during the chain formation the con- 
ductance ranges from 0.6 Gq to 1.1 Gq, exhibiting a long 
flat plateau close to 1 Gq. With respect to the conduc- 
tance value, this is in good agreement with experiment^. 
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FIG. 2: (Color online). Formation of Au chain in the [100] 
direction at T = 4.2 K. The upper panel show the strain 
forces as a function of the elongation of the contact. In the 
lower panel we show the conductance (solid line), the radius 
of the minimum cross-section (dashed line), and the channel 
transmissions (lines in the bottom of the graph). The vertical 
lines define regions with different number of open channels 
ranging from 8 to 1. Below the graph we show snapshots of 
the breaking of the contact. 



We observe that the conductance is mainly dominated by 
a single conduction channel. As in the case of a one-atom 
contact, the fact of having a single dominant channel is 
a consequence of the fact that Au is a monovalent metal. 
At the end of the stretching process, when the chain is 
made up of four or more atoms, we typically observe the 
appearance of a second channel and sometimes even a 
third one (see Fig. (21). Our analysis shows that these 
additional channels are due to the contribution of the 
d orbitals. With respect to the fluctuations of the con- 
ductance during the formation of the chains, we want to 
stress that they are not an artifact of our conductance 
calculations, but they are clearly related to fluctuations 
in the contact geometry. In Fig. (21 one can see that the 
structure in the conductance during the chain formation 
is accompanied by abrupt changes in the force values, 
i.e. by plastic deformations, which correspond to the in- 
corporation of another atom to the chain. We have cal- 
culated the conductance for all the configurations every 
5.6 ps, including configurations that could be unstable. 
In the experiment the stretching velocity is several orders 
of magnitude smaller and such unstable configurations 
are averaged out. 
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FIG. 3: (Color online). Snapshots of the chain formation 
of Fig. |5| At an initial time 0.0 ps the atom number 1 has 
three nearest neighbors and the atom 2 has four neighbors. 
At 1.4 ps the bond between atom 2 and atom 3 breaks. At 
25.2 ps the atom 2 is only connected to two atoms of the 
electrode. At 26.6 ps another bond breaks and atom 4 is now 
in a similar environment as atom 2 at the beginning. 
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FIG. 4: (Color online). Measured total conductance of a 
Au sample (top panel) and single-channel contributions (bot- 
tom panel) as a function of the electrode distance™" (opening 
curve). The vertical lines define regions with different num- 
bers of channels ranging from 6 to 1. The temperature was 
below 100 mK. For experimental details see Refs. |4a.l4^. 



IV. COMPARISON WITH THE 
EXPERIMENTAL RESULTS 

In this section we show experimental results obtained 
following the technique of Refs. p^ l42.l43l |. where part 
of the data has been already presented. In these exper- 
iments it is possible to extract the full information of 
the individual transmission coefficients. This is done in- 
ducing proximity superconductivity in a Au contact and 
analyzing the super-conducting current-voltage charac- 
teristics^^'^^.. We refer the reader to Ref. [43 for further 
details. Fig. ^ shows the conductance and the transmis- 
sion coefficients of a Au contact as a function of the elon- 
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FIG. 5: (Color online). Measured total conductance of a 
Au sample (top panel) and single-channel contributions (bot- 
tom panel) as a function of the electrode distance*" (opening 
curve). The vertical lines define regions with difi^crent num- 
bers of channels ranging from 8 to 1. The temperature was 
below 100 mK. For experimental details see Refs. |4a.l4^. 



gation at 100 mK. We can see how the channels disappear 
one by one as the contact evolves. Notice that there is 
a very long plateau of about 1 nm length, which corre- 
sponds to the formation of a chain. In this plateau the 
conductance changes between 0.3 Go and 0.9 Go and it 
is dominated by a single channel. Notice also the abrupt 
changes of the conductance, which are most likely due to 
incorporation of new atoms to the chain. 

In Fig. |S1 we show another experimental conductance 
curve. As it can be seen, there is a long plateau (about 
0.7nm), where the conductance is close to 1 Go. Naively, 
one would expect the conductance to be dominated by 
a single conduction channel. However, the analysis re- 
veals the presence of two channels, and only in the very 
last stages it reduces to one channel. In the two channel 
regime of the experiment in Fig.jSjthe second channel has 
a transmission mainly below 0.2 Go and the total conduc- 
tance is slightly below 1 Go. In our simulations we typ- 
ically observe two channels in the last stages for geome- 
tries in which there are two atoms in the narrowest part 
displaced with respect to each other (see Fig. EJa)). If 
there is only a single atom, the second channel gives only 
a small contribution to the conductance (see Fig.|H|^b)). 
As soon as we observe the dimer (see Fig.|B{c)), the con- 
ductance is largely dominated by a single channel. Our 
analysis of the character of the channels suggests that for 
a single central atom like in Fig. Etb) the second chan- 
nel is due to direct tunneling between the electrodes. In 
the dimer configuration the distance between the elec- 
trodes is considerably larger, which leads to the disap- 
pearance of the second channel. The current then flows 
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FIG. 6: (Color online). Transmission of the second channel 
during the last stages of the breaking of a typical An contact 
in the [100] direction (T=4.2K). The transmission of the first 
channel is about 0.7 Go. 



only through the s orbitals of the central atoms. 

The connection between the atomic structure and the 
number of channels described above is supported by the 
experimental curve of Fig. taken during the closing of 
a contact. In contrast to Fig. 0] and Fig. and all the 
simulations, where the contacts were stretched (produc- 
ing a so called opening curve), in Fig. [7| the two pieces of 
a nanocontact were brought into contact again after dis- 
ruption (thereby producing a so called closing curve). It 
seems very unlikely to produce a dimer by reconnecting 
a teared nanocontact, since the atoms collapse back after 
rupture into the bulk. In fact, one expects a single atom 
coupled to the two pieces of the nanocontact as smallest 
contacts, i.e. one expects two channels at the beginning. 
Indeed, the typical closing curve (Fig.[7|) starts with two 
open channels. 



V. MINIMUM CROSS-SECTION HISTOGRAM 

In order to discuss in more detail the mechanical prop- 
erties and stability of the contacts we have calculated 
histograms of the radius of the minimum cross-section 
for the three main crystallographic directions of the fcc- 
lattice. The results are summarized in Fig. |H1 In this 
figure we have normalized the radius to the radius of 
an ideal linear chain Rq = 1.268 A with nearest neighbor 
distance (lattice constant/ V2). The most prominent fea- 
ture of the histograms is the presence of peaks, especially 
for the [100] direction. These peaks suggest the existence 
of exceptionally stable radii. This is relatively surprising 
at low temperatures (4.2K), where the gold atoms do 
not have enough kinetic energy to explore many config- 
urations in order to minimize the surface energy. 

It is important to remark that these peaks generally 
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FIG. 7: (Color online). Measured total conductance of a Au 
sample (top panel) and single-channel contributions (bottom 
panel) as a function of the electrode distance^ (closing curve). 
The vertical lines define regions with different numbers of 
channels ranging from 2 to 5. The temperature was below 
100 mK. For experimental details see Refs. [43.14^. 



do not appear at the values ^/n ■ Rq (n — 1,2,3, . . .), 
and therefore we cannot infer that the cross-sections as- 
sociated with these peaks correspond to multiple of the 
cross-section of a single-atom contact. Furthermore in 
the three cases the peak close to 1 is extremely pro- 
nounced. This is due to the fact that usually in the 
last stages a dimer configuration or a chain of atoms is 
formed. 

The next obvious question is whether this peak struc- 
ture is an evidence of shell effects as observed for al- 
kali metalsii and suggested more recently for gold con- 
tacts22*24. For shell effects the peaks are expected to be 
equidistant when plotted as a function of the radius R of 
the wire^^. To test this idea in the inset of Fig.|Hl we plot 
kpR as a, function of the peak number for the histogram 
of the [100] direction. Here, kp is the Fermi momen- 
tum, which is assumed to take the free-electron value of 
1.21 X 10* cm~^. As one can see, the peak positions in- 
deed follow a straight line with slope 0.60 ± 0.02. This 
value of the slope lies in between the values observed in 
Ref. |2J| for electronic shells (1-02) and for atomic shells 
(0.40). From our analysis we cannot draw any conclu- 
sion on the appearance of shell effects in our simulations 
for several reasons. First, the choice of the peak posi- 
tions is not absolutely unambiguous, and second, a simi- 
lar analysis for the histograms in the directions [110] and 
[111] does not show such a clear linear relation between 
the radius and the peak position. It is worth stressing 
that even though the peak structure of these histograms 
cannot be easily interpreted in terms of shell effects, it 
clearly shows that particularly stable configurations do 
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FIG. 8: (Color online). Histograms of the radius of the mini- 
mum cross-section R, which is normalized in units of the ra- 
dius of an ideal linear chain of gold atoms _Ro ~ 1.268 A. The 
temperature is 4.2K. The different panels correspond to differ- 
ent crystallographic orientations. The number of simulations 
used to construct the histograms is 50 for each directions—. 
The inset in the upper panel shows the peak positions in the 
[100] histogram, converted to kpR, as a function of the peak 
index. The arrows in the upper panel indicate the peak posi- 
tions used to construct this plot. The last peak of the [100] 
histogram is artificially generated by the starting configura- 
tions. 



exist, which most likely correspond to meta-stable ge- 
ometries which are formed at low temperatures. 



VI. CONDUCTANCE HISTOGRAM 

In this section we want to address the main issue of this 
work, namely the origin of the peaks in the conductance 
histograms. Are these peaks simply due to the peaks 
in the minimum cross-section histograms as suggested 
in Ref. [2^? In order to answer this question, we have 
collected the conductance calculations^ for all our MD 
simulations in the [100] direction at 4.2K. The results are 
displayed in the bottom panel of Fig. |5| As one can see, 
the conductance histogram has a pronounced peak close 
to 1 Go and two further maxima above 2 Go and 3 Go. 
As we show below, the peak close to 1 Go is highly corre- 
lated with the first peak in the corresponding minimum 
cross-section histogram of Fig. |31(upper panel), and it is 
therefore a consequence of formation of dimer configura- 
tions or chains of atoms. However, the other peaks are 
not reflected in the conductance histogram in contrast 
to the postulated correlation between these two types of 
histograms for Al22i. We explain this fact as follows. As 
we pointed out in section III, the minimum cross-section 
is not the only ingredient that controls the conductance, 
and the geometry of the narrowest part and the disorder 
in the contact play also an important role. This means in 
practice that contacts with different radii can have simi- 
lar values of the conductance. In particular, this is true 
for small contacts like the one we have analyzed in this 
work. Of course, the situation could be different for larger 
contacts, where semi-classical arguments are believed to 
provide a good description. 

Let us now illustrate the correlation between the con- 
ductance histogram and the minimum cross-section his- 
togram. For this purpose, we show in Fig. |3 both 
histograms in a color/grey scale representation. The 
color/grey scale of the counts shows how many counts 
of a certain conductance value corresponds to a certain 
region of the radius histogram. Notice that the peak in 
the conductance histogram close to 1 Go is clearly due 
to the contribution of contacts with the cross-section of 
an atom, i.e. single-atom contacts or chains of atoms. 
The other minimum cross-section regions overlap in the 
conductance histogram. 

Another important piece of information can be ob- 
tained by analyzing the individual transmission coeffi- 
cients. Fig. ^1 shows the mean valued of the channel 
transmissions < T„ > as a function of the total conduc- 
tance for the conductance histogram of Fig. |3| Notice 
that below 1 Go the conductance is clearly dominated by 
a single channel. In particular, for a total conductance of 
1 Go the second channel gives only a contribution about 
0.036 Gq. This is due to the fact that this region corre- 
sponds mainly to dimer configurations or chain of atoms, 
where all the current proceeds through the 6s orbital of 
gold, as explained in section III. 

In the two channel regime of the experiment in Fig. [S] 
the second channel has a transmission mainly below 
0.2 Go and the total conductance is slightly below 1 Gq. 
Similar values are obtained in the simulations. For ex- 
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FIG. 9: (Color online). The lower panel shows the conduc- 
tance histogram for An in [100] direction for T = 4.2 K. We 
focus on the region G < 5Gq after the initial relaxation of the 
ideal contacts—. The color/grey scales indicate how many 
counts of a certain conductance value correspond to a certain 
region of the histogram of the radius of the minimum cross- 
section. This latter histogram is shown in the upper panel 
with the different regions defined*^ with different color/grey 
scales. In order to make it easier to distinguish the different 
colors in grey scaled print the colors are numbered. 



ample, in Fig. ^] we observe at a total mean conduc- 
tance of about 1.1 Go, i.e. in the two channels regime 
(the third channel has only a mean transmission of under 
0.03 Go) J a mean transmission of the second channel of 
about 0.18 Go. For higher conductance values, the chan- 
nels are not completely open, which explains the absence 
of conductance quantization in the histogram. 



VII. CONCLUSIONS 

We have presented a comprehensive theoretical study 
of the mechanical and electrical properties of Au nano- 
contacts at low temperatures and provide a comparison 
with experimental results on the transport channels. Our 
main goal was to analyze how the interplay of these two 
type of properties is reflected in the conductance his- 




0.5 



1.5 2 2.5 

conductance (G^) 



3.5 



FIG. 10: (Color online). Mean value of the transmission coef- 
ficients (Tn) (n = 1 until n = 6) as a function of the conduc- 
tance in units of Go for Au in [100] direction for T = 4.2 K, 
i.e. corresponding to the histogram of Fig. |5] The error bars 
indicate the standard deviations of the numerical results (if 
enough data points are available) . 



tograms. We have shown that the histograms of the mini- 
mum cross-section show a peak structure, which suggests 
that there are contact configurations which are particu- 
larly stable. Semi-classical arguments, which are mainly 
based on the Sharvin formula or modifications of it, in- 
dicate that the conductance is only controlled by the 
ratio between the radius of the contact and the Fermi 
wave length. With these arguments one would expect to 
see also a series of peaks in the conductance histograms, 
which are correlated with the peaks of the histogram of 
the minimum cross-section. However, in our calculations 
this is not the case, and only the peak close to 1 Go is 
clearly correlated with the first peak of the histogram of 
the minimum cross-section. In other words, we find that 
(i) the minimum cross-section is not the only ingredi- 
ent that determines the conductance and (ii) the peak at 
1 Go is due to the formation of single-atom contacts and 
chains of atoms. Our theoretical conductance histogram 
is in qualitative agreement with the experimental findings 
at low temperatures. We note however, that we cannot 
expect more than qualitative agreements due to the dif- 
ferent stretching process realization in the experiments 
and in the simulation. 

We have also found the following interesting results. 
In our simulations the most typical geometry in the last 
stages of the breaking of the contact is a dimer config- 
uration where two Au atoms are facing each other (see 
Fig. ^ . Our analysis of the conduction channels shows 
that the conductance is dominated by a single channel 
in the last plateau, which corresponds to a single-atom 
contact or a chain of atoms. Moreover, the channels dis- 
appear one by one as the section of the contact is reduced. 
Both observations are due to the fact that Au is a mono- 
valent metal, where every atom contributes with a single 
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valence orbital to the transport. These finding are in 
good agreement with the experimental observations (see 
Fig. 1415(1 . Furthermore the simulations indicate a small 
contribution of the second channel to the total conduc- 
tance if only one single atom is coupled to the electrodes. 
This finding can explain the number of channels of the 
opening curve of the experiment (see Fig.0). In addition, 
as shown in Fig. [S] and in agreement with our simula- 
tions, the last conductance plateau in the experiment is 
not always carried by only one single channel, as is often 
assumed. 

Of course, it would be desirable to extend our analysis 
to higher temperatures to investigate the appearance of 
shell efFects2i24. On the other hand, it would be interest- 
ing to compare the results for Au with results for other 
materials, in particular multivalent ones, to explain the 
differences in the conductance histograms. Both type of 
studies are currently under progress. 
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